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Abstract 

In this note we establish the existence of ground states for atoms within several 
restricted Hartree-Fock theories. It is shown, for example, that there exists a 
ground state for closed shell atoms with N electrons and nuclear charge Z > N —1. 
This has to be compared with the general Hartree-Fock theory where the existence 
of a minimizcr is known for Z > N — 1 only. 

1 Introduction 

Computations of the electronic structure of atoms and molecules in quantum chemistry 
in general rely on numerical solutions of simplified versions of the quantum many-body 
problem at hand. Among those, the Hartree-Fock approximation often serves as a 
starting point for more accurate approximations such as multi-configuration methods, 
see for example [8l[T6]. In the simplest version of Hartree-Fock theory the energy is 
minimized with respect to antisymmetric tensor products of orthonormal one-electron 
orbitals, the so-called single Slater determinants, and further restrictions are imposed 
in numerical procedures implementing this variational problem [4j. In any case the 
question arises whether a minimizer exists. This paper is concerned with several re- 
stricted Hartree-Fock theories for atoms where the one-electron orbitals are products 
of space and spin wave functions. For each of the considered restrictions we investigate 
the existence of a minimizer both for neutral atoms and positive ions, as well as for 
simply charged negative ions. 

The existence of a minimizer in the general Hartree-Fock ( GHF) theory for neu- 
tral atoms or positive ions was first established in 1977 by Lieb and Simon [11]. No 
constraints were imposed in their work besides the orthonormality of the one-electron 
orbitals. In the meantime there has been remarkable further progress in the study 
of the variational problem for the Hartree-Fock energy functional. It is known, for 
example, that there exists a sequence of critical points for this functional p^, and con- 
vergence properties of various algorithms used for the approximation of critical points 
were investigated in [U EJ [9] . 

The main concern of this article is the minimization of the Hartree-Fock energy 
functional under additional constraints. Our general assumption is that the one- 
electron states are products of space and spin functions. First, we treat the restricted 
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Hartree-Fock (RHF) functional for closed shell atoms with prescribed angular momen- 
tum quantum numbers. Second, we drop the latter requirement, i.e. we consider atoms 
with an even number of electrons, where only pairs of spin up and spin down electrons 
with the same spatial function occur. The corresponding energy functional will be 
called spin-restricted Hartree-Fock (SRHF) functional. It turns out that there exists a 
ground state in both cases, if Z > A^ — 1, where Z denotes the nuclear charge and N the 
number of electrons. The existence of a ground state in the case Z = N — 1 reminds of 
the well-known stability of closed shell configurations in chemistry. Third, we look at 
another restricted Hartree-Fock functional, which is called unrestricted Hartree-Fock 
(UHF) functional in the chemical literature, and must not be confused with the GHF 
functional. In the UHF setting, the spatial functions corresponding to spin up resp. spin 
down functions are chosen independently from each other, but still are assumed to have 
prescribed angular momenta. In this case a ground state exists if Z > — 1, and we 
provide sufficient conditions under which this is also true if Z = — 1. For example, 
there exists a ground state for Z = N — \ m. the spinless case (i.e. if all spins point in 
the same direction) with two angular momentum shells ii =0, H.2 > 0. 

For certain closed shell atoms (e.g. He, Ne) it is known that the minimization 
problems for the general and restricted Hartree-Fock functionals coincide, if Z ^ N 
On the other hand there are also cases where they differ [14] , see [7] for an explanation 
of this fact. Nevertheless, the restricted ground states are always critical points of the 
GHF functional. This is due to the fact that the considered constraints do not require 
additional Lagrange multipliers in the Euler-Lagrange equations. Thus, this paper also 
establishes the existence of critical points for the GHF functional in the case Z = N — 1. 
To our knowledge, the only previous result providing the existence of critical points for 
the GHF functional in the case Z = — 1 is given in the paper [5] of Cances and Le 
Bris, which in fact even holds for arbitrary Z > 0. But in general, the critical points 
constructed in their paper only correspond either to local (not global) minima or saddle 
points. 

In the literature the existence of minimizers for restricted Hartree-Fock functionals 
has previously been studied for special cases. Based on Reeken's paper [13] on the 
solutions of the Hartree equation, Bazley and Seydel [2] proved the existence of a 
minimizer for the spin-restricted Hartree-Fock functional of Helium (A^ = 2), which is 
given by the restricted Hartree functional. For this functional it is known that there 
exists a minimizer even if Z = 1 = N—1, see [12^ Theorem II. 2]. In our paper we extend 
this result to arbitrary numbers of filled shells. Lieb and Simon generalize their GHF 
existence result [TT] to certain restricted situations in [lOj, but their theorem does not 
cover the restrictions discussed in this paper. Though, this article has been strongly 
inspired by their work [11]. In [12], Lions treats restricted Hartree-Fock equations, 
which arise as the Euler-Lagrange equations of the RHF functional. He proves the 
existence of a sequence of solutions to these equations provided Z > N. Lions' proof 
relies, however, on the unproven assertion that all eigenvalues of a radial Fock operator 
are simple. His approach is motivated by the paper of Wolkowisky [T7] who shows the 
existence of solutions for a system of restricted Hartree-type equations. A numerical 
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approach to restricted Hartree-Fock theory may be found in the book of Froese Fischer 
[6j. Finally, we mention the article of Solovej [15], where he proves the existence of a 
universal constant Q > so that there is no GHF minimizer for Z < N — Q. This 
establishes the ionization conjecture within the Hartree-Fock theory. The question 
whether or not there is a GHF minimizer for Z = N — 1 is open. 

The paper is organized as follows: In Section [2] we introduce the restricted Hartree- 
Fock functional for closed shell atoms with prescribed angular momentum quantum 
numbers and prove an existence theorem for minimizers of this functional. The Sec- 
tion [3] is devoted to generalizations of the RHF existence theorem to the SRHF and 
UHF functionals. A derivation of the RHF functional in the closed shell case can be 
found in Section [H Finally, there is an appendix containing technical lemmas. 

Acknowledgment. The author thanks M. Griesemer for drawing his attention to the 
problem and for helpful discussions. The author is supported by the Studienstiftung 
des Deutschen Volkes. 



2 Minimizers for Closed Shell Atoms 

The simplest Hartree-Fock approximation for atoms consists in restricting the admis- 
sible A^-electron states to the set of single Slater determinants, which are of the form 

A • • • A V7iv)(xi, . . . ,XAr) = — L= ^ Sgn(cr)99^(i)(xi)...V9^(Ar)(j:iv), (1) 

where Sn denotes the symmetric group of degree A^, sgn((7) is the sign of a permutation 
cr, and ipi, . . . , ipj\[ denote orthonormal L^(M^; C^)-functions. It is well-known, that the 
energy of an atom with nuclear charge Z and electrons in the state ([1]) is given by 
the general Hartree-Fock ( GHF) functional 

e^'i^.- . . , = E / |V^/ - ^1,./ + 1 // a. dy (2) 

where 

N N 

T{x,y) := ^ipj{x)'^j{y), p{x) := ^ I'/'jl^)!^ 

i=i i=i 
denote the density matrix and the electronic density, respectively. We write x = 
(x, /i) G M'^ X {±1}, J dx refers to integration with respect to the product of Lebesgue 
and counting measure, and |x — y| = |x — y|. 

Given a closed shell atom with sq G N shells of prescribed angular momentum 
quantum numbers ^i, . . . ,iso £ Nq, we impose the following form on the one-electron 
orbitals 

¥'jm<7(x,/u) = ^^^I^^Y£^.„(x)5^^, j = 1,. . . ,so, m = . . cr = ±1, (3) 

where the radial functions fj are in L^(M_|_) and 

(/^, fj) ■■= I Tifj dr = 6,j, if = £j, (4) 
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in order to ensure the orthonormality of the functions Here Yim denote the usual 
spherical harmonics. The Hartree-Fock energy of the Slater determinant built by 
the orbitals ^ is given by the restricted Hartree-Fock (RHF) functional (compare 
Section HI: 

e^'^'ih, . . . , /.o) = 2f:(2^, + 1) ( / l/jl^ + - f I// dr) 

2^.^i \JJ{R+r max{r,s} 



-2fj{r)Ms)Ue^e,ir,s)Mr)fjis)drdsj. (5) 

The integral kernels Ui-i^, appearing in the last term on the right-hand side are given 
in (j25p . We shall only need their properties collected in Lemma |5.1[ 

Let ifg(M+) denote the completion of Cg°(M+) with respect to the ff-'^(]R-|_)-norm. 
The RHF functional ([5]) is bounded below, if the functions /i, • • • , /sg are in ifQ(M+) 
and obey the constraints see Lemma 15.21 We define the RHF ground state energy 
by 

EiN, Z) = inf{f . . . , /.Jl/i, . . . , fs, G i^oHK+), (/., /,) = % if = ^,}, (6) 

where the dependence of E{N,Z) on ii,. . . ,isQ is suppressed. The main question of 
this paper is whether the infimum in ^ is actually a minimum. 

If there exist minimizing functions fi,...,fsQ obeying the constraints of ([6]), then 
they are solutions of the corresponding Euler-Lagrange equations, which we may as- 
sume to have the form (see Remark (b) below) 

Hiji = eifi, i = l,...,so, (7) 

with radial Fock operators given by 

He, = -52 + M^i±il_^ + 2C/-K4, i = l,...,so, where 

{Uf){r) = f2{2e, + l)[ M^dsfir), 
^ 7lR+ maxjr, sj 

{Kef){r) = V(2£, + l)/,(r) / j;{s)f{s)Uu^{r,s)ds. 
j=i 

We suppress the dependence of the operators U, Ki and thus Hi- on the functions 
/i ) • • • 5 /so • The Euler-Lagrange equations ([7]) , called Hartree-Fock equations, form a 
set of So coupled non- linear eigenvalue equations for the functions /i , . . . , /^p . 
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Remarks. 

(a) By Lemma 15. 2( the operators Hi. are symmetric semi-bounded operators on 
C^(M+). Therefore, minimizing functions fi, . . . , fs^ obeying the constraints 
of dl]) are in the domain D{Hi.) of the Friedrichs extension of -fT^,., which is 
contained in i?Q(M-|_). 

(b) The Euler-Lagrange equations for minimizing functions /i , . . . , /^y obeying the 
constraints of ([6]) are given by Hi-fi = Ylj^ijfj^ where the sum runs over all 
indices j with £j = £j. Since the functional £^^^ is invariant under unitary 
transformations of the subspaces of (M+ ) spanned by all radial functions fj with 
equal angular momentum quantum numbers, the minimizing functions /i, . . . , /^p 
can always be chosen as eigenfunctions of the radial Fock operators. 

(c) The constraints of ([6]) may be relaxed without lowering the ground state energy, 
more precisely E{N, Z) = E{N, Z) for 

E{N,Z) = mf{£^^^{fu...Jso)\fi,---,fso(^H^{^+),{fiJj)=Oife,=ej 

and i ^ j, \\fi\\ < 1 for all i}. (8) 

This can be seen using similar arguments as for the general Hartree-Fock func- 
tional in \12\ section II. 2]. The following theorem shows that the relaxed mini- 
mization problem always possesses a minimizer. 

Theorem 2.1. Let sq £ N, ii, ■ ■ ■ ,iso ^ Nq, and Z > 0. Then, there exist functions 
/i, . . . , /so G //o(^^+); which minimize the RHF functional ([5]) under the constraints 

{fi, fj) = Oifii = £j and i 7^ j, 
II /ill < 1 for alii. 

Moreover, fi £ D{HeJ, Heji = Sifi, and: 

(i) Either £i < or fi = 0. £i < implies \\fi\\ = 1. 

(ii) If Z > N - 2{2ii + 1), then fi / 0. 

If Z > N — 1, then \\fi\\ = 1 for all i = 1, . . . , sq. 

(Hi) If Z > N — 1, then £i < and \\fi\\ = 1 for all i = 1, . . . , sq. 

Remarks. 

(a) Theorem 12. II (ii) shows that for Z = N — 1 there always exists a normalized mini- 
mizer 

{oj.£RHF_ In this case we do not know whether or not < 0. Nevertheless, 
it is clear that Z > iV-2(2^j + l) always implies E{N, Z) < E^-^^N -2{2£i + l), Z) 
for alH = 1, . . . , So, where E^^\N — 2{2£i + 1), Z) denotes the minimal energy in 
the case where all electrons of the i-th. shell are dropped. 
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(b) In general Hartree-Fock theory it is known that the minimizing functions can 
be chosen as eigenfunctions to the lowest eigenvalues of the corresponding 
Fock operator. Moreover, there is a gap between the occupied and unoccupied 
eigenvalues [Ij. It would be interesting to know whether similar results hold also 
in the restricted Hartree-Fock theory, where, unfortunately, the method of [Ij is 
not applicable. 

Before turning to the proof of Theorem 12.11 we introduce the following notations 
that will be used throughout the paper: 

r> := max{r, s}, r< := min{r, s}, for r, s > 0. 

We write £^^^{fi, . . . , fi, . . . , fs^) to denote the restricted Hartree-Fock functional 
where the electrons of the i-th shell are dropped. The following lemma exhibits the 
dependence of £^^^{fi, . . . , fi, . . . , fso) on /j, and will be crucial for the existence of a 
minimizer in the critical case Z = N — 1. It follows easily from the definition of 

if we set P,{r,s) := (21, + l){2r^' - Ue^e,{r,s)). 



Lemma 2.2. Let sq G N, 



more, let i G {1, . . . , sq} and let H^^^ denote the Fock operator where all electrons of 
the i-th shell are dropped. Then: 



G No, Z > and fi, , 



,fso G 



Further- 



£ 



RHF 



ifli • • • 1 fit • • • 1 /so) 



cHHf ( f f f ■ 

^ Ul) ■ ■ ■ ) Ji^ ■ ■ ■ 1 J So. 

+ m^ + l){fi(^fi\Pi\fi 



where Pi{r,s) = Pi{s,r) and 

2£i + l 
max{r, s} 

Furthermore, for all A > 0; 



< P^{r,s) < 



max{r, s} ' 



+ 2{2l, + I] 
^fi), 



r,s > 0. 



f^ 



fi 



(9) 



^^^^(/i, ...,/,,..., /,J + 4(2^, + 1)5 Re {h \Hi^\f-: 



fi + Sh 



ISO) 



+2i2ii + l)6^{h 



H. 



h 



\{fi\H,^ \ fi) + Re{h® h\P,\ fi® fi, 
+{fi®h + h® fi \Pi\ fr ® /i)) + 0((53) 



(10) 



for 5^{). 



Proof of Theorem \2.1l First, we give a proof of the existence of a minimizer for the 
relaxed minimization problem, which proceeds the same way as in the paper of Lieb and 
Simon [TI]. £^^^{gi, . . . , ^r^g) is bounded below independently of gi, . . . , gso £ -ffg (ffi+) 
with \\gi\\ < 1, see Lemma 15.21 Thus, let g^\...,gi^^ be a minimizing sequence for 
the relaxed minimization problem ([8|). Again by Lemma 15.21 (^fj"^)^^^ is bounded 
in //q(M+). Hence, there exist weakly-f^Q (1R+) convergent subsequences g^"''^ 
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n — )• oo), j = 1, . . . , So- Fix i G {1, . . . , sq}- Without loss of generality we may assume 
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that gi, . . . ,gk^ are all functions Qj with 



i. The matrix M := {{gj, gk))j,k=i,...,k, is 



hermitian and obeys < M < 1 (c.f. [Ill Lemma 2.2]), so there exists a unitary ki x ki 
matrix U with the property U*MU = D, where D is a diagonal matrix with eigenvalues 
in [0, 1]. If we define fj = Yl'k^i Ukj9k, j = '^,---,ki, then {fj, fk) = Xj6jk, < Xj < 1. 
It is easy to see that £^^^ is invariant under such transformations. Thus, transforming 
each subspace of functions with equal angular momentum quantum numbers in this 
way, we obtain functions /i, . . . , /sp with {fi, fj) = 0, if ii = ij, i j, \\fi\\ < 1 for all 
i. Furthermore, fi, . . . , fs^ minimize £^^^ , because 

EiN,Z) < ^:«^^(/i,...,/,J=f^^^(5i,...,<7so) 



< 



lim inf £ 

ra— )-cxD 



1 iJso ) 



E{N,Z), 



; fs, 



where we used Lemma 15.21 By further transformations we can achieve that /i, 
are eigenfunctions of the operators H£.. 

(i) Let fi and assume that > 0. Then, by (fTOl) with A = and h = fi, the 
energy decreases if we decrease the norm of fi. Let < and assume that ||/i|| < 1. 
Then, the energy is decreased by increasing the norm of fi. 

(ii) We prove the following more general statement: Let < /i < 1 and let 
Z > iV - 1 - (1 - fi){Aii + 1), then fi < \\fif < 1. 

There is nothing to prove in the case fi = 0. Therefore, let > and assume that 
< n. We show that there exists h G Hq(R+) with h A- fj, if Ij = ii, such that 

^''"'^{fl, ...J^ + 5h,..., fs,) < ...,fi,...Jso) 



for small S ^ 0, which contradicts the minimization property of /i, . . . , fsg. The depen- 
dence of the left-hand side on h G Hq{W^) is given by (jlOp with A = 0. The factor of 5 
in (jlOp vanishes since fi, . . . , fs, is a minimizer. Therefore, we only have to show that 
there exist infinitely many normalized functions h G Hq{E.^) with disjoint supports, 
such that the factor of 5"^ in (IIOD 



(^h l^lf I + {fi®h \P,\ f^®h) + if, h\P,\ h®fi)+ Re ih h\P,\ f, ® fi) (11) 

is negative. We may drop the Re-term because it becomes non-positive upon a suitable 
choice of the phase of h. Let J E C^(M+), supp(J) C [1,2], ||J|| = 1. Furthermore, 
we define Ji?(r) := R-^/^J{r/R) for R > 0, then supp(Ji?) C [R,2R], \\Jr\\ = 1, 
Jr e C^(M+). Using U{r) < r"! EjLi(2^j + 1) and Ki>0 (LemmaE^]), we see that 



Jr 



H 



Jr 



< 



Jr 



-92 + 



+ 1) 



Z N 
— + — 
r 



2i2ii + 1) 



Jr 



:i2) 



This inequality combined with the estimate for Pi in Lemma 
([TT]) with the choice h = Jr 



allows us to estimate 



{fi^JR\P^\f^i 
{fi®jR\Pt\jR 



Jr 
Jr) 



< 



< 



mf£ 



R 



'J i J 
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for i? — )• (X). The sum of the three terms on the right-hand side becomes negative for 
R large enough, because < /U, by assumption. This proves (ii). 

(iii) It suffices to show that ej < 0, j = l,...,so, see (i) and (ii). Assume that 
£i = 0. We show that there exists h G //q(M+), \\h\\ = 1, h ± fj, if ii = ij, so that 

pRHFff fi + f \^cRHF(f f f \ 

^ KJl-i ■ ■ ■ 1 r- ^) ■ ■ ■ 7 J SI) ) ^ ^ \J 1: ■ ■ ■ : J iy ■ ■ ■ : J so ) 

V 1 + 0^ 

for small 6^0. Again, the dependence on h of the left-hand side is given by (jlOp with 
A = 1. Since £i = 0, it suffices to show that the factor of (5^, which is the same as in 
(fTTj) . can be made negative by suitable choices of h. This can be done choosing the 
same scaled functions as in (ii), but now using Z > N — 1 instead of < □ 



3 Other Restricted Hartree— Fock Functionals 



Theorem 12.11 can be readily generalized to other restricted Hartree-Fock functionals. 
In this section we present analogous results for a so-called UHF functional as well as 
for a spin-restricted Hartree-Fock functional. In the former case we continue assuming 
that the electrons are in product states of space and spin but we drop the condition 
that the spatial wave functions for spin up resp. spin down electrons are equal in each 
shell with fixed angular momentum quantum numbers. More precisely, we consider 
electrons that are in states of the form 

/;(|xi). 



(/?jm-f-(x,/i) : 

where Sq , Sq G Nq, 



|x| 
//(|x|) 



y,.^(x)5^,+i, j = i,...,s^, m = -e'^, +q 



-1 ' ■ 



G No, and for all v G {«,/?}, i,j G {1, . 



(13) 
(14) 



/rGi7i(K+), (/r,/n = 5,„ if^r = A- 



The corresponding Hartree-Fock functional, which is called unrestricted Hartree-Fock 
(UHF) functional, takes the form 



£ 



UHF/ fa 



f fa fa . fP fP \ 

Wl 1 • • • ) /sgi /l > • • • 5 / jsj 



f/3 



ue{a,l3}j=l 



. rji'' + 1) 



■'J 



*0 *0 



iJ&{ot,l3} 3,k=l jr=l fc=l 

Here we use the shorthand notations 



(15) 



{2t^ +l){2t^ + 1)1 f^®f^ 



3 k 
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Given u G {a, /?} and ^ € No we introduce Fock operators 

HI := -dl+i{i + l)r-^ - Zr~^ + U - K"^, where 

for / e L^(R+). Again these operators depend on the functions /f , . . . , /^^. Using the 
same methods as in the proof of Theorem 12. H the following existence theorem can be 
proved: 

Theorem 3.1. Let s^,s^ e No, £f, . . . ,£^^,1^, . . . j'^f, G No, and Z > 0. Then, there 

exist functions /f , . . . , /f , . . . , /^^ G -ffg (R+), which minimize the UHF functional 
(fT5]) under the constraints: for all v G {a, /?} and j G {1, . . . , Sq} 

(/r,/;) = ifei = t), i^j, 

ll/ril < I for alii. 
Moreover, f- G D[H-.), H^.f^^ = e'^ f- . 

(i) Either e\ <0 or fj^ = 0. < implies \\f[\\ = 1. 

(ii) IfZ>N- {2£^ + 1), then / 0. 
IfZ>N-l and / 0, i/ien \\f^^\\ = 1. 

(Hi) If Z > N -I, then e\ < and \\f^\\ = 1 for all v G {a,/3}, i = 1,. . . 

iiemarfe. 

(a) We do not know, except for the case where £ = 0, whether the occupied eigen- 
values of the corresponding Fock operator are the lowest eigenvalues or whether 
there is a gap between occupied and unoccupied eigenvalues. 

(b) In general, Theorem 13.11 does not imply the existence of UHF minimizers in the 
case of Z = A'^ — 1. Nevertheless, in the special case where all spins point in the 
same direction (i.e. the spinless case) the following existence result holds true. 

Corollary 3.2. Let s% G N, = 0, and let £f = 0, l^, . . . > with 

If Z = ^j=2(^^f + -'-) '^'^^ N = Z + 1, then the UHF functional (jl5p has a minimizer 
under the constraints (/°, /°) = 5ij for all j = 1, . . . , Sg with ii = £j. 
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Remark. The condition of Corollary 13.21 always holds in the case of two shells Sg = 2, 
£° = 0, £^ > 0. 

Proof of Corollary \ 3.^ Theorem 13.11 yields the existence of /f,...,/^ G //q(IR+), 
which minimize (fTSl) under the constraints (/f , /") = if = and ^ / j, || < 1 
for ah i. Clearly, ||/2°|| = • • • = || = 1 by (h). Observe that 

^^^^(/r,...,/.p< inf ^^^^(<7,0,...,0) = -^, (16) 



Il9ll<l 



and on the other hand 



gUHF^Qj-a^ ja^ > f£« + 1)2 ~ T 1 ^ ^ ( F+T ) 



> 4- 

where we dropped the electron-electron energy and estimated the remaining terms by 
the hydrogen ground state energies in the first inequality, and used the condition on 
Sq in the second inequality. Assume, that (/f l^^o l/f ) ~ 0, then 



, ...,f%) = ^^^^(0, f^, f%), 



because ^^^^(/f , ■ ■ ■ , f^^) = £^"^{0, /a", • • • , /"^) + {f^\H^\f^), which contradicts 
(USD and (HZ]). Therefore," (/f |i^o l/f ) = ll/f f °< 0, which implies < and thus 
ll/f II = 1- □ 

Another frequently used model for atoms with an even number of electrons is the 
spin-restricted Hartree-Fock (SRHF) model [4]. It emerges from the RHF model in 
Section [2] by dropping the prescribed angular momentum quantum numbers. More 
precisely, for an atom with atomic number Z and N = 2n we impose the following 
form on the one-electron orbitals 

ipiai:>c,n) = (Piix)6a^,, z = l,...,n, a = ±1, 

where ifi G H^CR.^) and cpj) := jj^g ^ipj dx. = 6ij. Then the restricted Hartree-Fock 
functional reads 

£"""'(Vi,...,V„)=2X; I |Vv,(x)|2- 

2JJ |x-y| |x-y| 
Here the electronic density matrix and the electronic density are given by 

n n 

■^(x,y) = ^fi{x)ifi{y), p(x) = ^|v9i(x)|2. 

i=l i=l 
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The corresponding Fock operator is given by 

m J |x-y| 

where {Kip){x.) := J "^^^^^yi^^ dy. Using similar ideas as in the proof of Theorem l2. li the 
following existence theorem holds true for the spin-restricted Hartree-Fock functional: 

Theorem 3.3. Let Z > and N = 2n. Then, there exist functions (pi, . . . ,ipn G 
H^(M.^), which minimize the SRHF functional (fT8|) under the constraints 

{ipi,ipj) = Oifij^j, 
\\(Pi\\ < 1 for all i. 

Moreover, ipi e D{H) = H'^{R^), Hipi = Si^pi, and: 

(i) Either £{ < or ipi = 0. £i < implies \\(pi\\ = 1. 

(a) If Z > N — 2, then ipi ^ for all i = 1, . . . ,n. 
If Z > N — 1, then \\(pi\\ = 1 for all i = 1 . . . ,n. 

(Hi) If Z > N — 1, then £i < and \\<fii\\ = 1 for all i = 1, . . . ,n. 

Remark. For this spin-restricted Hartree~Fock functional the minimizer exists for all 
Z > N — 1. Again we do not know whether or not £j are the n lowest eigenvalues of 
H, although there seem to be no numerical counterexamples [Ij. 



4 Derivation of the Closed Shell Energy Functional 

For the reader's convenience we give here a self-contained derivation of the restricted 
Hartree-Fock functional ([5]). For this purpose, we begin with a lemma that will be 
useful for the calculation of the electron-electron interaction energy. 

Let Pi denote the ^-th Legendre polynomial. We remark that for x, y G and 
£ G No the addition theorem 

2£ + 1 

Yumuf) = -^p^^ ■ y) (19) 

m=—£ 

holds, where x • y is the usual scalar product of two vectors in W^. 

Proposition 4.1. Let i,L e No and M G Z, \M\ < L. Then for all r,s > and 
X G §2. 



1 f m-y)YLM{y) , ^ ^ (L i nV minjr,.}" 

4^ U |rx - .y| ^"^^^ = ^^^^("^ 1^0 oj max{r, .j'^^^ ' 



2 
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Remark. An easy consequence of this proposition is that for all i, i' € No 

1 f m-y)PA^-y) y^' (l i' kV min{r,5}^ 



This is seen by multiplying ([20|) with iLAf(x), integrating over with respect to x 
and summing over M = —L, . . . ,L. 

Proof. Assume first that r ^ s. For fixed x G the series expansion 

1 ^ — 



rx — sy\ r> \ / 



converges pointwise for all y G and thus in ^^(S^) because ^^=0 ( ) • y) is 

bounded uniformly in N and y. We get 

Pei^ • y) 



rx — sy, 

' ' n=0 



1 / \ " 



> n=0 ^ ^ fc=^-n V / 



where we used the addition theorem 

2 

I Pk{z). 



i+n / 

Pn{z)Pi{z) = ^ (2A; + 1) 

k=\l-n\ \ 



kin 




The addition theorem ()19p allows us to compute 

1 / PH-yji-.MK),^,.) 



47r y§2 |rx — sy 

00 / \ n f-+n 



n=0 \ \ / m=—k 



00 



yLM{^)Y. 



L £ n\ min{r, s}' 



n=0 



0/ max{r, s} 



71+1 ■ 



The desired equation for r ^ s now follows from the fact that the Wigner 3j-symbols 
vanish unless |-L — £| < n < L + i. The case r = s can be derived from the above 
result by choosing a sequence r„ 1 s. Clearly, r-^—^ t .,-,1 for all y € \ |x) and 
j -.^-. is integrable with respect to y G S^. Hence Lebesgue's Dominated Convergence 
Theorem may be used to see that the formula is also true for r = s. □ 

Let us turn to the derivation of £^^^ . If /x , . . . , /^^ are in H^iM.+), then the 
functions tpjma defined by (l3|) are orthonormal in L'^ (W^ ; C"^) , and ipj-ma 
by Hardy's inequality 

^^^dr < 4 / |/'(r)pdr (22) 
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for / E iJQ(R+). Using the addition theorem (fT9]) . the corresponding density matrix r 
and electronic density p take the form 



so 



r(x,y) = 5^.M,2^^- rn r^P^,-(x-y), (23) 



47r X 



Here we abbreviate x := x/|x| for all 7^ x G M'^. If the general Hartree-Fock 
functional ([5]) is evaluated at the functions (/'jmo-! the only term, which is not trivially 
computed, is the exchange term: 

12 ■50 



\^-y\ j^^i i^^y Ak+p 

r ,^ P,^.(x-y)P4(x-y) 
X / do- x,y — ^— — . 

y(s2)2 FX - sy\ 

Using (|21|) . the form of (jS]) follows from the choice 

f fcV min{r,4fe 

5 Appendix 

Lemma 5.1. Let 1,1' £ No, andr,s > 0. T/ien the functions Uu/ defined by (|25p o6e?/; 
("f/ij Uu>{r,s) = Ui'i{r,s) = Uui{s,r), 
(U2) < Uu,{r,s) < max{r,s}~i. 



(U4) For all g G ffg(M+) the integral kernels g{r)Ui£i{r, s)g{s) define non-negative 



Hilbert-Schmidt operators on 



Proof. (Ul) and (U3) are obvious from the explicit representation of Ui£'{r,s) and 



\ / 2£ + l 

(U2) The positivity of Uw is clear, the upper bound can be proved using ([2T]) . ([19] 
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and Cauchy-Schwarz: 

(4vr)^ J(S2)2 |rx - sy| 



m=—t m' - 



(2£ + l)(2f + \y(s2)2 |rx-sy 



(§2)2 |rx-sy| 



(47r)2 7(g2)2 |rx — sy| ' max{r, s} ' 

where we used 2ab < + b"^ and (|19|) in the last inequahty. 

(U4) The integral kernels K{r,s) := g{r)Uiii{r, s)g{s) are in L^(M^) by (U2) and 
by Hardy's inequality (j22p . which shows that the corresponding integral operators are 
Hilbert-Schmidt. Moreover, let 

93m(x,/u) := ^4^^£m(x)5;,,+i, m = -£,... 
|x| 

and 



( \ \- ( \ X X 2^ + l g(|x|) ff(|y|) 

Given / G L^(M+), we define 



I ^ ={21+1)11 W)K{r, s) /(.) dr d.. 



then 



The last equality can be computed using (pOj) and (j25|) . Hence, the non-negativity of 
the integral operator corresponding to K follows from the non-negativity of the term 
on the left-hand side. □ 



Lemma 5.2. (i) For all f G H^{R+) and e > 0; (/, i/) < e\\f'\f + 



2 , 1 II ^||2 



(ii) Let So G N, ^i, . . . ,4o G No, Z > 0, /i, . . . , fso G H^(^+), and e > 0. Then: 



f«^^(/i,...,/,J>2^(2£, + l) 



so r rv 

'l|2 ^ II ^ ||2 



(l-Ze)||/j„ 



(^mj Let So G N, £i, . . . ,4o G No, and /i, . . . , /sq G i?o(M+). Then for all I G No-' 

so 

0<ir,<[/< J](24 + l)(||/^f + 
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(iv) Let E Nq. Then the following maps are weakly sequentially continuous on 
H^{R+) resp. H^{R+) x H^{M+): 

if, 9) ^ (/®5'|max{r,s}^^| 
if, 9) ^ {f ®9\Uu'\g® f)- 

(v) The functional £^^^ is weakly sequentially lower semicontinuous on xf^-^H^^Rj^) . 

Proof, (i) and (iii) follow easily from the Cauchy-Schwarz and the Hardy inequalities 
(f22D . (Ul), (U2), and (U4). To prove (ii) fix j, A: E {1, . . . , sq]. Using Cauchy-Schwarz, 
(Ul) and (U2) we obtain 



< 



fj{r)fk{s)Ui^i,^{r,s)fk{r)fj{s)dr ds 

\f,{r)\^\fk{s)\''U,^,,{r,s)drds 
\fo{r)\^ms)\^Ut^t,{r.s)drds 



\fk{r)\''\fj{s)\^U,^,,{r,s)drds 



< 



\Mr)\^\fk{s)\' 
max{r, s} 



dr ds. 



so 



Therefore, 

i=i 

The claim now follows immediately from (i). 

(iv) Let fn^f weakly in ffQ(R_(_). Due to the Rellich-Kondrashov theorem, 
converges to / uniformly in M+. To prove the weak continuity of the Coulomb potential 
we first use 

K/n,i/„>-(/,i/>| < K/n-/,i/n>| + |(i/,/n-/>| =(*) + (**). 

For R > we obtain using Cauchy-Schwarz and Hardy's inequality (j22p 



(*) < 







^ l/n(r)-/(r)||/„(r)| ^ 1 
ar+ ^ 



< 



r 

R \ 1/2 

|/„(r)-/(r)|2dr 



\fnir)-fir)\\fnir)\dr 



R 







dr 



1/2 



+ ^ 1 1 /n / 1 1 1 1 /n I 



< 2VR sup {|/„(r)-/(r)|}||/;|| + -(||/„|| + 11/11) 



re{0,R) 



R 



Since ||/;;||, 11/11, and ||/„ 



are uniformly bounded in n, we can first choose R large to 
make the second term small, then choose n large to make the first term small. (**) can 
be estimated analogously. The weak continuity of the other maps can be seen with a 
similar decomposition argument as shown above for the Coulomb potential. 



(v) Let ff-^ fj weakly in H^{R+) for j = 1, . . . , TV. Clearly, 



fj 



_Q2 ^ inil 



/,) <liminf(/ 



,(n) 



2 , 



An) 
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since fj fj in iJQ(]R+) implies drfj drfj in L^(M_|_) for the first term, and 
using the lemma of Fatou for the second term. The remaining terms of S^^^ are 
weakly sequentially continuous as shown in (iv). □ 
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